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Abstract 

We suggest a class of direct sample estimates for the two-point quadratic statistical functions 
of multidimensional data, which includes: estimates of the sample autocovariance function (AcF), 
sample mean square increment (also, structure) function, and the estimate of the power spectrum. 
The central estimate for the class is the sample AcF, which is constructed as to represent the finite 
Fourier transform of the periodogram estimate of the spectrum and is positive semidefinite. The 
estimate explicitly account for the anisotropy of the fields in all spatial directions and is illustrated 
on two examples: the morphology of the Grab nebula and the surface roughness generated as a 
solution of the anisotropic Kardar-Parisi-Zhang equation. We also provide an expression of the 
covariance of the sample AcF in the case of data assumed to be drawn from a two-dimensional 
Gaussian field with a known mean. 
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I. INTRODUCTION 



Despite the increased availabihty of multidimensional data, it appears that not much 
attention has been devoted to the problem of adequate and accurate direct estimation of 
simple quantities such as two-point, quadratic statistical functions. The latter include au- 
tocovariance (AcF) and autocorrelation functions, and the mean square increment (also 
structure) function. Among the monographs we are aware of, a good exception is ^1] , where 
an estimate of the two-dimensional AcF is discussed in some detail. Specifically, for a given 
Ni X N2 matrix of real data fxiX2J ^1 = 1;2, ...A?"!, X2 = 1,2, ... N2, the autocovariance 
function is estimated in volume 2, chapter 9 of |1] by 

^ Nj—ui N2—U2 

= ^ ^ ^ ] ^ ] fxi,X2 fxi+Ul,X2+U2^ (1) 
^ ^ Xl=l X2 = l 

where u = (mi,U2), fxiX2 = [1x^x2- Pj with / denoting the sample mean, / = 
{I/N1N2) T.xl=i T.xl=i f-m The expression is vahd for < Mi < A^i-1 and < M2 < A'2-l 
and is extended to the third quadrant using the AcF property of being even function. For- 
mally, the latter replaces ui and U2 in (PJ) by their absolute values. Note that this leaves 
A^{u) undefined in the second and the fourth quadrants. We shall call this estimate "the 
standard" 2d sample autocovariance function (SAcF) estimate. 

For a homogeneous (stationary) random 2d field, f{x), the spectral representation theo- 
rem asserts that its AcF and spectral density function (or simply the spectrum) are related 
by the Fourier transform, 

(1 \ 2 +00 +00 
' xi=—oo X2=~00 

where k = {ki,k2). We use calligraphic letters to distinguish the "true" spectra and AcFs 
from their estimates. The inverse transform of Q reads: 

/•TT l*TT 

A{x) = / S(k)e-''^-^cPk. (3) 



Equations Q and Q apply to a random field, which depends on discretely valued vector 
variable. It is desirable to have estimates that are also related by a (discrete and finite) 
Fourier transform. The simplest direct estimate of the spectrum is provided by the peri- 
odogram 

^ Afi-l A^a-l 

^^^^ ^ (27r)2A^iA^2 ^ ^ fxiX2fyiy2e "^^'^ (4) 

2:1, yi = l X2, 1/2 = 1 



The periodogram, referred to as "raw" (unsmoothed) estimate, is a basis for variety of more 
sophisticated spectral estimates. Hence, in accordance with we would like to have an 
AcF direct estimate, which for all < ImiI < A^^i — 1 and < \u2\ < N2 — 1 satisfies 



A{u) = I{k)e-'''-''d'k. (5) 

J —TV J —TV 

The standard estimate, A^,, satisfies (0), in the domains where it is defined. However, it is 
easy to see that if the validity of is extended to the second and the fourth quadrants, 
Eq. (0) does not hold. The latter follows also from the developments presented in the next 
section. 

Another property of the theoretical AcF that must be shared by an estimate is the 
property of positive semidefiniteness. Let arbitrary non-zero {Ni x A^2) matrix 

of real numbers. Then the AcF of a 2d stationary field satisfies 

ax^,X2ayuy2-^{xi-yuX2-y2)>0- (6) 

xi,yi=l X2,s/2=1 

It easy to see that A* does not satisfy © in the second and the fourth quadrants. We stress 
that the inequality © does not bear theoretical importance only; it ensures the positiveness 
of spectral estimates based on an AcF estimate (indirect estimates) P], see also Eq. Q 
below. 

In this paper we obtain a 2d AcF direct estimate, denoted hereinafter by A{u), which 
satisfies (0) in all four quadrants, see the next section - Eq. (jHJ. The estimate is identical to 
A*('u) in the first and the third quadrants but differs in the second and the fourth quadrants. 
The estimate A{u) is positive semidefinite and leads to new estimates of both the mean 
square increment (structure) function and the power spectrum. We briefly discuss these 
estimates also in the second section and provide generalization of (jH)) for arbitrary dimension. 

One of the most important advantages of the new estimate lays in the fact that it is ca- 
pable to capture the anisotropy of the data in arbitrary spatial direction. The latter dictates 
the choice of the illustrations and applications we consider here; yet another application — 
a study of an YBCO thin film morphology using (jSI) — can be found in P]. In the first 
of the two applications, we study the AcF of images representing the morphology of the 
Crab nebula. The random field for these images is the light intensity in a recorded pixel. 
The anisotropy in this case is determined either by the direction of the expansion of the 



supernova ejecta or by the interaction of synchrotron nebula with the ejecta 



m. 



In the third section we compute the covariance of A{u) for dimension d = 2 and under 
the simphfying assumptions that the observations are drawn from a Gaussian random field 
and are already adjusted to have zero mean. Our second application involves the anisotropic 
Kardar-Parisi- Zhang (AKPZ) equation j^Q,!^, considered in section four in somewhat more 
detail. The equation pertains to growth of vicinal surfaces and the anisotropy arises from 
the different rates of growth along and across the average steps direction [6]. We study 
how this anisotropy imprints on the shape of the AcF by numerically solving the AKPZ 
equation on a large lattice and then taking a smaller size images rotated on various angles 
with respect to the axis of anisotropy. A summary of our results and main conclusions are 
presented in the last section. 



II. SAMPLE ESTIMATES OF THE MULTIDIMENSIONAL ACF 

In order to obtain a direct estimate of the 2d sample AcF that corresponds to I{k), 
we simply need to rearrange (jlj) to a 2d discrete Fourier transform. We begin by chang- 
ing yi = Xi — Ui and reversing the order of Xi and Ui summation. This breaks up 
the sums with respect to xi,yi to two double sums: Yl!ll=i'Yliyl=i = Yl!ll=iYTu^l-^i = 
Elli=-(7Vi-i) ES"" + T.ul=lT.xl=i+u,- Next, we shift the summation xi ^ Xi + Ui in 
the second term only obtaining the intermediate result: 



^ N2 ( Q Afi-|Mi| 

(2nVN^N', 2^ \ 2^ 2^ Jx^,x2Jx^+\u^iy2^ 

y ' ^ ^ X2,y2=l I ui=-(Ni~-l) xi=l 

Ni-1 Ni~\ui\ 

1 \^ f f p-i[kiui+k2{x2-y2)] 

^ / J / J JX\^-\u\\,X2JXx,y2 ^ 

2:1=1 

In analogous manner we deal with the X2, and 2/2 summations obtaining four terms, which 
then can be combined into two pairs arriving at the following 2d discrete Fourier represen- 
tation of the periodogram 

^ A^i-l N2-I 

^(^^) = 7^ E E me-^'- (7) 

^ ' ui=-Ni+lu2=-N2+l 
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In ((7j), we introduced the function 

Ni-\ui\ N2-\u2\ 



A{u) 



Jf]^ •^^i'^^ /xi + |«i|,X2+|«2| for m ■ M2 > 

^ Xl=l X2 = l 

^ A^i — N2 — \u2\ 

Y Y /^l-^2 + |«2| fxi + \^i\,X2 for Ml • M2 < 

^ ^ Zl=l X2 = l 



which represents the new estimate of the sample AcF in dimension two. It differs from the 
standard estimate (Q) in the second and fourth quadrants only. We stress that by the virtue 
of its construction, it is ^ but not the standard estimate that is related to the periodogram 
by©. 

A different way to obtain A{u) is to substitute I{k) directly in Eq. ^ and perform the 
integration with respect to k. Then to use the obtained pair of Kronecker symbols to carry 
out two of the summations carefully accounting for the limits of summation. The latter 
depends on the position of u, which produces (jHl). 

The estimate A{u) is positive semidefinite function. This follows from the way it was 
obtained but it is instructive to show that the positive semidefiniteness virtually forces the 
form of the SAcF, Eq. (jH)). To demonstrate this we use a 2d generalization of a proof of 
semidefiniteness due to McLeod and Jimenez, see chapter 6. Let exi,x2 be a 2d white noise 
with zero mean and variance I/N1N2. Form a the field gxi,x2 = fk,i£xi-k,x2-u here 

the centered values of the sample, fk^u are considered as coefficients of a moving average- 
like 2d field. Since gxi,x2 is stationary, its theoretical AcF must be positive semidefinite. 
A straightforward calculation shows that this function is given by (jH)). One last remark 
regarding (jH)): it easy to see that the estimate A{u) can be computed numerically using 
the FFT-based algorithm and the codes given in [9], chapters 12 and 13; however, by first 
extending the data twice in both dimensions and assuming that fxx,x2 = ^/ either Ni < 
Xi < 2Ni, or N2 < X2 < 2N2 (double length zero padding). 

Estimate (jH)) can be used as a basis for estimates of both sample power spectrum and 
sample mean square increment function. One can obtain an estimate of the spectrum from 
(jH} by employing it in a 2d generalization of Grenander and Rosenblatt formula [lfll |. 

Afi-l N2-I 

^(^) = 7^ E E Hu)Aiu)e-^'-^, (9) 

ui=-Ni+lu2 = -N2+l 
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where the function w{u) is termed "lag window". The general properties and specific exam- 
ples of w{u) can be found in we note here that obviously it must be an even function 
w{—u) = w{u). The statistics of this estimate will be presented elsewhere, note however that 
as in the Id case, (0) provides a smoothed compared to the "raw" periodogram estimate. 



The mean square increment (structure) function of /(x), B{u), is defined jll|, |3| by 
B{u) = E [{f{x + u) — /(5))^], where E denotes the ensemble averaging, and is related 



to the AcF by I3{u) = 2 {A{Q) — Mu)) , e.g. '13]; for an estimate that corresponds to the 



3,13. 



An estimate corresponding to ^ should be modified 



standard estimate of A{u) see 
in the second and the fourth quadrants as: 

^ Ni-\ui\ N2~\u2\ 

B{U) = {fx^,X2+\u2\ - fx, + \nilx2T , for Ui-U2< 0. (10) 

^ ^ Xl = l X2 = l 

In addition of being a positive function, we remark that B(u) does not involve the sample 

_nn . . 

mean of /(x) and thus is free from a source of bias brought up by / [2, |16|. Modifications 
analogous to (fTUj) are apparently due for the estimates of the generalized structure functions 



used to infer multifractal scaling, see Refs. 
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Estimate (jH)) can readily be generalized to arbitrary dimension d. To shorten the nota- 
tions, consider a multi-index with d components p = {pi,pi, . . .pd) each taking a value of 
either or 1. Let = indicates x^, whereas pk = 1 indicates Xk + \uk\, k = 1,2, .. .d, and 
let p designate the multi-index whose components are all different from the components of p; 
e.g. if p = (Oil), then p = (100). Then the (i-dimensional SAcF estimate can be expressed 
as 

Xl = l Xd = l 

where p = (0, 0, . . . , 0) pertains to the first and the 2"^^^ + 1-th hyperquadrants, p = 
(0, 0, ... , 1) to the second and 2°^"^ + 2-th hyperquadrants, and so on. In general, the 
SAcF in hyperquadrants k and 2'^^^ + A; is expressed by ()11|) with a multi-index p, which is 
the binary representation of number k — 1. 

To illustrate the difference between the standard and the estimate (jHI), we provide plots 
of both estimates for pair of images representing regions of the Crab nebula. The images 
were selected from a color image taken from the Hubble Space Telescope (HST) The 
color image was created as a weighted sum of three narrowband filters centered at 5012 
A, 6306 A and 6732 A and comprises 24 individual Wide Field and Planetary Camera 2 



exposures ^. We converted the color image into gray scale (in the range of 0-256) image, 
i.e. fx,y in this case is proportional to the light intensity per pixel and represents the 
morphology of the nebula; we added a circle indicating the position of the Crab pulsar. 
The first of the selected images, shown in Fig. Wip) top panel, is located close to outer 
rim of the nebula dominated by the expanding ejecta jj, |5|. The coordinates of its right 
bottom corner are, RA: 02000 =5:34:40.5 and Dec: ^2000 =21:59:39.1. The image extends 
(46.2x49.7) arcseconds corresponding to {Ni = 453) x (A^2 = 487) pixels. By just inspecting 
the image the anisotropy is not easily recognizable, however, due to the outward expansion 
of the supernova remnant an anisotropy roughly across the radial direction (direction to the 
pulsar) should present in the morphology. 

The second image has coordinates: ^2000 =5:34:29.0, ^2000 = 21:59:10.0 and extends 
51.8x51.8 arcseconds, (504 x 504) pixels — Fig. ^b), bottom panel. It is from a region 
where the synchrotron nebula (upper left sector of the image) interacts with the denser 
ejecta creating "filaments" . The latter are attributed to a magnetic Rayleigh-Taylor (R-T) 
instabilities The major axis of anisotropy in this region should, in general, be expected 
along rather then across the direction to the pulsar. 

An overall linear (planar) trend, g = s^a + ss6 with origins of coordinate systems at the 
left-upper corners of the images, is removed before both SAcF estimates were computed. 
Carrying out this procedure is important since the trend by itself produces anisotropy in 
the SAcF. In the case of images presented in Fig. ^^b) the linear trend is rather small; 
the estimated slopes (in units grayscale/pixel) are sdec = —0.027, sra = 0.0071, and 
■sdec = 0.101, Sra = 0.031, for first and second images, respectively. The standard and the 
estimate (jH)) for the first image are presented as gray scale plots with superimposed level 
contours in Fig. ^c) and Fig. ^d), respectively. For the sake of completeness of the plot the 
standard AcF is extended to the second and the fourth quadrants, hence the specific "rose" 
appearance of the AcF. The anisotropy of the supernova remnant structures in this part 
of the nebula is clearly recognized from the plot of A{u) with the major axis of anisotropy 
having an angle of about ~ —25° with respect to the horizontal axis. This angle should be 
interpreted as the average front of the local expansion, refer to Fig. ^a). Another quantity 
that characterizes the anisotropy is the aspect ratio, 7, the ratio between the characteristic 
sizes of the nebula structures along the minor and majors axes of anisotropy. The latter 
sizes are defined by the correlation lengths of SAcF in the respective directions. We evaluate 
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FIG. 1: Standard AcF estimate and estimate Q for two regions of the Crab nebula. Panel (a) 

shows the location of the regions; the circle indicates the position of the pulsar. Panel (b) shows 

the two regions zoomed by an identical factor. Panels (c) and (d) show the standard AcF estimate 

and the estimate respectively for the upper image in (b). Panels (e) and (f) - the same the 

bottom image in (b). The morphology principle axes of anisotropy are drawn in plates (d) and (f). 
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these lengths crudely by assuming, somewhat arbitrary, that the Id principal profiles of AcF 
are represented by random processes with finite domain (band-limited) spectra. U sing the 
expressions for the correlation length obtained for this class of random processes in j20|, |21|, 
we infer 7 ^ 0.57. 

For the second region, the major axis of the anisotropy is at angle of -(9 ~ —59.5° with 
respect to the horizontal axis, which as should be expected is roughly in direction to the 
pulsar, refer to Fig. IH^a) and Fig. [Iff). The aspect ratio in this case is 7 ~ 0.875. 



III. COVARIANCES OF THE 2D SAMPLE AUTOCOVARIANCE FUNCTION 

In this section we obtain expression for covariances of 2d SAcF — the estimate (jHI), 
evaluated at two points u = (mi, U2) and v = (f 1, ^2): 

cov [A{u),A{v)] = E[A{u)A{v)] - E[A{u)]E[A{v)]. (12) 

The covariance has both theoretical as well as practical importance for determining the 
confidence intervals in the AcF estimate. The expression will be derived under the simpli- 
fying assumptions that fxi,x2 is a Gaussian random field with zero mean (or that the mean 
is known and subtracted). It is immediately seen that irrespective to which quadrant u 
belongs, 

m^)] = (1 - ^) (1 - ^) ^("). (13) 

where as before A{u) denotes the true autocovariance function of f{x). Eq. (fTT^ is identical 
to the ensemble average of the standard estimate and shows that (jHI) is a biased estimate. 
The bias, however, is small for large samples especially for small \u\. We turn now to the first 
term in ()12j). Reckoning with A{—u) = A{u) and the symmetry under the exchange u ^ v, 
we need to consider three different cases only: (i) n G I quadrant, -u G I quadrant; (ii) u G 
II quadrant, t? G II quadrant; and (iii) n G I quadrant, -u G II quadrant. The calculations in 
all three cases are closely similar; below we illustrate them for the case (ii) only. Inserting 
the pertinent for this case AcF expressions from (jH)) and using that for a Gaussian field the 
four-point function can be expressed as combinations of products of two two-point AcFs we 
have 
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E[A{u)A{t 



^ Ni-\ui\ N2~-\u2\ Ni-\vi\ N2-\v2\ 
^ 2 xi=l X2=l yi=l s/2=l 

+A{yi-xi,y2-X2 + \v2\-\u2\)A{yi-xi + \vi\-\ui\,y2-X2) (14) 
+A{yi-xi + \vi\,y2-X2-\u2\)A{yi-xi-\ui\,y2-X2 + \v2\ 

The first of these terms does not depends of the summation indexes and cancels out with 
the second term in (fT^ exactly, refer to ()1H|1 . In the remaining two terms we perform the 
indicated change p = y — x and again reverse the order of the summation. This allows to 
carry out the summations with respect to both Xi and X2 explicitly, noting in the process 
that we need to distinguish the case ImiI > |fi| from the case > The result is a 

product of two trapezium shaped functions, which involve two parameters, c and d, 



W{c,d-p) 



0, 

1 - c + p, 
1 — c — d, 
1 - c - p, 
0, 



p< -(1- 
[l - c) <p< -d 
-d <p< d 
d <p< 1 — c 
1 — c <p , 



c) 



(15) 



see also Fig. |21 The parameters c and d are subject to the conditions Q<d<l — c<l. 
Next, introducing a := {u + v)/2\ b := {u — v)/2 and shifting the summation indexes 



Pi and p2 simultaneously according to p = (gi 



Ml 



Ni,q2 + U2 — N2), we arrive at the 



following expression for the covariance valid when both u and v are in the second quadrant: 



cov 



A{u),A{v) 



2Afi+2ai-l 2N2-2a2-l 



qi = l 52 = 1 

A{f+ a) A{f- a) + A{f+ b) A{f- b) 





M 


ri ^ 






r2 \ 




iVi 


iViy 


N2 


N2 



(16) 



where r = (gi — Oi — A^i, ^2 + «2 — ^2) has been introduced. 

Similar expressions are obtained in cases (i) and (iii) above. Finally, if we define 0,1^2 '■= 
(1^1,2! + |?^i,2|)/2 and 61^2 := (|wi,2| — |'yi,2|)/2, all three cases can be combined into the 
following compact form: 
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FIG. 2: Graph of the window function W{c,d;p), Eq. H15|) . used in the expression of the sample 
AcF covariances. 
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ai \bi\ Ti 
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gi=l 92=1 
A{r + a) A{f- a) + A{f + b) A{f - b) 



(17) 



with the general definition r = (gi + ai — A^i, ga + ^2 — A^2)- This expression is vaUd for 
arbitrary positions of vectors u and v. 

The important for the practice variances of the SAcF, can be obtained from the expres- 
sions of the covariance. Setting v = u'ui (fTTjl we have 



var [A{u) 



1 



Af2-|M2|-1 



1 - 



l^il + bi 

iVi 



1 



\U2\ + \P2\ 



X 



pi=-Ni-\\ui\H p2=-N^U2\n 
•^{Pl,P2) + A{pi + \Ui\,P2 + \u2\) A{pi-\Ui\,P2-\u2\) 



(18) 



Note that we went back from summation with respect to qi and q2 to summation with 
respect to pi and p2, which results in symmetric about zero limits of summation. 
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IV. APPLICATION TO THE ANISOTROPIC KPZ EQUATION 



The anisotropic Kardar-Parisi- Zhang (AKPZ) equation has been introduced in an at- 
tempt to model the growth on a vicinal substrates jy]. Adatoms that migrate towards the 
steps and attach to them, have lower probability to desorb compared to those migrating 
parallel to the steps. This effectively induces different rates of growth along and across the 
steps and violates the rotational symmetry of the KPZ growth process j^j. Let h = h{x,t) 
be the height of the growing surface at point x = {x, y) and time t. If one chooses the 
x-coordinate along the direction of the steps, the AKPZ equation takes the form: 

dth = v^dlh + Vydlh + ^ {d^hf + ^ {dyhf + r/, (19) 

see also 0]. In this equation: and Vy are coefficients of the curvature terms associated 
with desorption, A^; and \y are coefficients related to growth rates normal the surface, and 
Tj = ri{x, y, t) is a Gaussian white noise, E[ri{x, y, t)rj{x' , y' , t')] = 2D6{x — x')6{y — y')6(t — t'). 
The equation has been studied by D. E. Wolf using one- loop, renormalization-group (RG) 
approximation Some of the obtained results have later been confirmed by numerical 
simulations [2^. To recap what will be needed here, let Tj, = Uy/u^ and rx = Xy/Xx and let 
both Xx and be positive. In this case the AKPZ surface grows with an exponent identical 
to the surface generated by the isotropic KPZ equation, referred to as algebraic roughness. 
As the morphology evolves, the nonlinear parameters A^^ and Xy, as in the case of isotropic 
KPZ are not renormalized, whereas and Uy take effective values such that r^, = rx (a fixed 
point to the djTiamical renormalization flow equations) . This means that in this case the 
anisotropy of the surface is of the simplest kind - elliptical anisotropy - and therefore might 
be taken as a benchmark for testing statistical methods characterizing anisotropy. 

The numerical simulations were carried out using the Amar and Family numerical 
a Q. b.oad,. .peaM„. .elude. .e.ea.,n. of .Ke e,ua.o„ and e.- 

ploying the standard discretization for the derivatives. Two comments are in order. First, 



in contrast to 



23|, we choose rescaling that leaves the equation manifestly anisotropic: 



X 1-^ ^/l^x, y (-^ \fv^y-, h (2i'nc/Xx)h, and t] t-^ y^Dlvjc]. Second, the discrete analog 
does not adequately represent the continuous AKPZ equation j^], however, since more ac- 
curate difference scheme are not known for dimensions higher than one 0, , we employ 
here the standard discretization. What is more important within the scope of this study. 
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FIG. 3: Gray scale image of a morphology obtained by numerically solving the anisotropic KPZ 
equation. The image represents the central (512x512) part from the entire (1024x1024) simulated 
surface. The conditions at which the simulation is carried out are discussed in the text. 

the discrete equation inherits the elliptical anisotropy of the original AKPZ equation. 

The elliptical anisotropy can be discerned even by visually inspecting of the simulated 
morphology, see Fig. El The picture represents AKPZ surface generated with parameters 
Vx = ^) = 0.3, = 10, \y = 3, and D = 0.2; (to skip the "transient" time for the system 
to come to the RG fixed point, we have chosen r,^ = rx = 0.3 at the outset). The simulation 
is carried out on a square lattice with side of L = 1024 and for T = 2 x 10^ time steps of 
At = 0.001. The surface height range is given in units of lattice spacing set to unity. 

In a typical experimental circumstances, the axes of anisotropy are rarely known and 
need to be inferred and quantified from a image of the morphology To reckon with 

this we "record" smaller, (512x512), images, which are rotated at angles ipo = 0°, 10°, 
30°, and 60° with respect to the x-axis of the simulated surface. The picture in Fig. El 
represents the image for tpQ = 0°. The images for ipQ ^ 0° are obtained using a simple, 
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FIG. 4: Sample AcF for an image of the AKPZ simulated morphology illustrated in Fig.|21 however 
rotated on angle 30°. 

based on the four nearest neighbor points interpolation. Then for every image we compute 
the AcF estimate (jH)), an example of which for ipQ = 30° is shown in Fig. 01 Taking a more 
systematic approach, rather than the correlation length, we consider sections of AcF defined 
by aA{0) < A{u) < {a + Aa)A{0) for several levels < a < 1 and a fixed width of Aa = 0.04. 
We project the values of A{u) within each section on the (ui, U2) plane and fit these points 
by an ellipse. The direction of the axis of asymmetry and the aspect ratio is evaluated from 
the parameters of these ellipses. In the actual fits we have used four levels of a: 0.2, 0.4, 
0.6, and 0.8. 

The obtained results show a discrepancy typically within ±2° from the expected direction 
of the anisotropy. In few cases only, all associated with the lowest level of a, the discrepancy 
is up to ±4°. More interesting are the inferred values of the aspect ratio. These are plotted in 
Fig.Elfor all four angles of rotation and AcF levels a and for two simulated surfaces, studied 
independently. The parameters of the first simulation are the same as those used to produce 



14 



o 

(0 



0.70- 

0.65- 

0.60- 

0.55- 

0.50- 



U 
0) 

Q. 0.45- 
(0 

0.40- 

0.35- 
0.30- 



□ 



X 

o 



X 

IT 



□ 


deg 


o 


10 deg 


X 


30 deg 


+ 


60 deg 



X 



0.548 



0.387 



— I ' 1 1 1 ' 1 ' 1 ' 1 ' 1 ' 1 

0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

AcF level a 



FIG. 5: Retrieved aspect ratio 7 from four AcF sections marked by levels a and for four angles of 
rotation of the recorded image, see the legends. Upper and bottom data represent results from two 
simulated surfaces, see the text. The horizontal lines indicate the expected values of the aspect 
ratio, 7 = ^Jvyjvx- 

Fig. ISl The lattice size, the time step, and total integration time for the second simulation 
are also the same, however, with: = 1, Uy = 0.15, A^; = 10, A^; = 1.4, and D = 0.2 
i^u = rx = 0.15). The values of 7 for the first simulation are grouped around 0.548, upper 
part of the figure, whereas for the second — around 0.387 bottom part of the figure. Both 
these values correspond to the respective ^/r^ values used in the simulations. To understand 
this, we rescale ()19|1 in a manner different from the one used prior to numerical integration; 
namely, x = y^z7~r, y = y/iA^y, h = {^^^J2D jvH'^vy'^ h, and r] = {^/2D/ux^'^i^l^^)rj, arriving 
at 

dji = dih + dfh + ^-^ {d£) ' + I {dyh) ' + ?j. (20) 

In the above equation, = j Vx'^ r\! ^ and Sy = {r\/ry) Sx and hence at the fixed point 

of the RG, Eq. (j20|l is equivalent to the isotopic KPZ equation. Therefore, the anisotropic 
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surface in this case is obtained by just rescaling the isotropic surface in x and y directions 



which as demonstrated in Fig. El is imprinted in the sample AcF, Eq. (jSJ. A somewhat 
larger discrepancy from ^Jr^ observed at the lowest level fits can be attributed to a greater 
relative variability of the SAcF. The latter can be estimated crudely from the variance 
fll8|) by substituting the SAcF for the unknown true AcF. For the two simulations used in 
Fig. El we obtain an increase from 12% at n = 0, up to about 23% for points u at which 
A{u) = 0.2A{0). As a final remark, the one-loop RG approximation of the AKPZ equation 
indicates that if Ay 7^ the surfaces in this class are characterized by two characteristic 
lengths, Ix and ly, even when the morphology has not yet evolved to the RG fixed point. 
In addition, the characteristic length scale linearly, ly ~ l^, Ref. [||. This suggests that the 
approach undertaken in this section may be suitable for characterization of more generic 
AKPZ morphologies. Further numerical simulations, however, are needed to confirm this 
assertion. 

V. CONCLUSIONS 

In this paper we have suggested an estimate for the autocovariance function (AcF) of a 
homogenous random field in arbitrary dimension d. The estimate, Eq. (jllll . is constructed 
as to represent the discrete and finite Fourier transform of the periodogram estimate of the 
field's power spectrum; it is identical to the standard AcF estimate in the first and the 
2d-i _j_ quadrants but differs in all other quadrants. As it should be, the estimate is 
positive semidefinite. On the basis of (fTT|) . we have suggested new estimates for the field's 
structure function and power spectrum. We have also derived expressions for the covariance, 
consequently for the variance of the AcF estimate in two dimensions under the simplifying 
assumption that the field is Gaussian and with a known mean. 

Perhaps the most important advantage of the new sample AcF over the standard estimate 
lays in the fact that it captures the anisotropy of the field in all spatial directions. The latter 
is demonstrated on two examples. The first involves the morphology of the Grab nebula 
observed by the Hubble space telescope. For sake of comparison we presented plots of 
the standard AcF as well. The second example involves surfaces simulated by numerically 
solving the anisotropic Kardar-Parisi- Zhang (AKPZ) equation and is considered in more 




respectively. The latter leads to aspect 
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detail. In particular, we have focused on the case Py/fx = \/^x, i-e. when the system is at 
a fixed point of the djTiamic renormalization group approximation for the AKPZ equation. 
In this case the characteristic lengths of the morphology are two and are determined by 
Ix = and ly = y/i^- Hence, the surface can be viewed as a simple benchmark for testing 
statistical methods that account for anisotropy. We have shown that one can retrieve both 
the direction and the aspect ratio of the anisotropy reasonably well from the estimate in 
two dimensions. This has been done on several sections of the AcF and on two independent 
realizations. 
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